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ABSTRACT 



This investigation considers a non-linear technique for 
extracting natural resonances from transient electromagnetic 
scattering responses of radar targets. These natural reso- 
nances represent the complex poles of the target's transfer 
function in the Laplace transform s-plane. The advantage of 
their use in identification is their dependence only upon the 
geometry and composition of the target and not upon the aspect 
and polarization of the incident signal. 

Based on recent theoretical research, a new- signal model 
has been developed for describing the form of the transient 
scattering response. This form precludes the optimal use of 
previous methods for natural resonance extraction based on 
Prony's algorithm. In this effort, a modified least-squares 
approach is taken which can accommodate the actual form of 
the transient signal. The algorithm is tested using simulated 
signals with various noise levels and conclusions are drawn 
concerning the viability of the method. 
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INTRODUCTION 



I . 

This investigation considers a non-linear technique for 
extracting natural resonances from transient electromagnetic 
scattering responses of radar targets. These natural reso- 
nances represent the complex poles of the target's transfer 
function in the Laplace transform s-plane. The advantage of 
their use in identification is their dependence only upon the 
geometry and composition of the target and not upon the aspect 
and polarization of the incident signal. 

Concepts and methodologies evolving from this idea have 
been developed based on the Singularity Expansion Method 
(SEM) by Baum [Ref. 1] in which the location of s-plane 
singularities characterizes a target. The direct extraction 
of poles (and residues) from the complex impulse response was 
initially attempted by Main and Moffat [Ref. 2]. These two 
early references spawned the idea of using natural resonances 
as the basis of non-cooperative target recognition (NCTR) . 

The application of the Householder orthogonalization 
method using a successive orthogonalization process and the 
eigenvalue method, that finds eigenvectors and eigenvalues, 
can determine the number and value of poles from the tradi- 
tional transient signal model of a scatterer if there is a 
sufficiently low level of noise in the data [Ref. 3]. But 
both methods show poor quality of performance if we apply 
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them to the new signal model, even were there is no noise 
assumed in the signal data. This new signal model, which 
represents the electromagnetic scattering transient response 
of actual targets, was recently derived by Morgan [Ref. 4]. 

As considered in [Ref. 2], the relationship between a 
scattering target and the singularities being extracted from 
the natural complex resonance signal has been emphasized in 
the context of the number of poles and the accuracy of poles 
to ensure the unique relationship between an object and its 
image. Namely, the accuracy of pole extraction plays an 
important part in setting baselines to check consistency 
between a received signal and a set of singularities. Also, 
[Ref. 2] describes the modeling problem of a signal such that 
we express a transient response signal in terms of complex 
exponentials and the way it can be manipulated by Prony ' s 
algorithm. This signal model (termed here as the traditional 
model) , was found not to represent the actual transient signal 
in the recent work by Davenport [Ref. 5]. In 1982, Manilha, 
in his work, commented that the new signal model might be 
partitioned into "early time signal" and "late time signal", 
each of which characterize the signal by the time instant at 
which the traditional model can be acceptable or not [Ref. 3]. 
This effort considers the development of a method that is 
capable of handling the "early time" whose time varying 
residues of exponentials do not become constant until the 
excitation of the incident field disappears. Such a signal 
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can not be adequately or accurately modeled by the tradi- 
tional model. 

It is the intention of this work to investigate a method 
that can handle the "early time signal" and show "Robustness" 
under relatively heavy noise pollution as well as improving 
the accuracy of the poles we are interested in evaluating. 

A non-linear parameter optimization approach, through the 
modified least-squares minimization method, has been tried. 

In particular, we have been interested in evaluating its 
performance in improving the accuracy of the poles extracted. 
The non-linear parameter optimization approach using itera- 
tions has contribued to improvement of accuracy and showed 
the possibility of application to the "early time signal" as 
was the final goal. We define the "new signal model" as a 
causal connection of the "early time signal" and "late time 
signal" throughout this paper. The new signal model will be 
presented in Chapter II, with both the description of the 
transient response mechanism and the way in which the signal 
model is constructed through the transient response mechanism. 
Chapter III briefly describes the problems associated with the 
traditional model, and with many of the current extraction 
methods in existence. Chapter IV is devoted to the descrip- 
tion of the non-linear parameter optimization through a 
modified least-squares method and the way in which the algo- 
rithm is implemented. Chapter V presents test results of an 
attempt to use this approach for pole extraction and its 
performance characteristics. Chapter VI contains a summary 
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of the simulation results, the problems encountered, the 
potential problems expected, and the possibility of 
application of this method. 

Computations were performed using the IBM 3033 system at 
the Naval Postgraduate School. The testing of concept feasi- 
bility was the major impetus behind this effort, with concerns 
for practical implementation in a "real-time" environment left 
for subsequent study. Considerations to the processing speed 
will not be addressed in this thesis. All the calculations 
were done in a double precision environment. 
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II. SCATTERING TRANSIENT RESPONSE SYSTEM 



A. TRANSIENT RESPONSE MECHANISM 

An electromagnetic wave indent upon a scattering body 
forces currents to be distributed such that Maxwell ' s equa- 
tions and the corresponding boundary conditions are satisfied. 
These induced currents produce a scattered field that can be 
defined at any coordinate in free space if the expressions for 
those currents are made correctly. Changes in the angle of 
incidence, the incident wave shape, the polarization mode, 
and the target geometry and composition affect the results of 
the evaluation of the scattered transient response . 

Basically, the singularities uniquely characterize a 
scatterer by their nximber and locations. The corresponding 
residues, which are interpreted as the weights of these poles, 
either can be viewed in the frequency domain or in the time 
domain . 

As shown in [Ref. 3], equation (2.1) represents the form 
of the scattered signal impulse response of the target, in 
digital form 

00 

H(k) = (kT) + I H. exp(kSjj^T) (2.1) 

^ m=l 

where H^(t) =0 for t > T0, T is the sampling interval and 
0 — 

the infinite sum represents the natural resonance response 

with complex S = o + in the left— half plane. The poles 

^ in m m 

constitute the parameters for NCTR of an object. The process 
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of extracting those parameters from the impulse response 
first collects the data points within a finite time period 
(time window) . 

The data collected is represented using equation (2.1) 
and equation (2.3) which partitions the signal into "early 
time" and "late time" components. 

00 M-1 

X(k) = I H(k) * I i(k) (2.2) 

k=0 k=0 

T0-1 M-1 oo M-1 

X(k) = I H(k) * I I(k) + I H(k) * I I(k) (2.3) 

k=0 k=0 k=T0 k=0 

early time data late time data 

where I(k) is the incident wave at t = kT 
T is sampling interval in time. 

We can write the equation (2.3) in the form of equation (2.4) 
or (2.5) . Both of these expressions display the essence of 
the new signal model which will be discussed in the next 
section . 

N 

X(k) = I h (k)exp(S^k) (2.4) 

n=l ^ 

where A (k) = for k > T0 
n n — 

N 

X(k) = E(k) + I exp(S^k) (2.5) 

n=l 

where E(k) = 0 for k ^ TJ? 

and A is constant for all k 
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T0 is the point at which the late time signal starts, so we 
consider that a time window might contain the data points 
ranging from the one before T0 to that after T0. If we 
perform processing on a data window containing the data 
points which were from that early time region, we may lose 
some of the poles of natural resonance what are likely to be 
seen in the late time region. The data points after K0 
presumably has all the information regarding poles, even 
though there are some extra poles that may be discriminated 
from the natural poles . These extra poles are introduced 
by the result of noise with relatively low SNR. Also, there 
is another risk to lose the natural poles having very high 
damping coefficients. There is no effective method yet found 
to compromise on having the point in which we can see all the 
poles with less extra poles under the higher SNR condition. 

Meanwhile, the discrimination of E{k) appears to us as 
the quantity to be handled to give improvement in accuracy 
of the natural poles. The next section describes the new 
signal model based on these aspects and the transient response 
mechanism. 

B. IMPULSE RESPONSE SIGNAL MODEL 

The late time impulse response can be expressed as a 
summation of exponentially damped sinusoids as in equation 
( 2 . 6 ) . 

N 

X (k) = ^ A exp(S k) , k=0,T,2T, . .. , (M-l)T (2.6) 

a n=l ^ 
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where T = At, and M is the number of sampling points. 

The signal model expressed in equation (2.6) is valid after 
the incident wave has completely illuminated a target and 
no forced response remains . 

As the mechanism of a transient response system is strictly 
governed by the physics and orientation of the target and the 
attributes of the incident wave, we can write the equation 
(2.7) as an implicit form of equation (2.2), assuming we 
strike the target with a plane wave impulse. 

°° T0-1 

X. (k) = I H(k) . I [u(k) - u(k-T0)] (2.7) 

° k=0 k=0 

Combining both equations (2.6) and (2.7), equation (2.9) would 
be sufficient if we assume N poles are present. 

X(k) = X^(k) + Xj^(k) 

00 T0-1 

X(k) = I H(k) • I [u(k) - u(k-T0)] 

k=0 k=0 

+ exp(S^k) , k=0,T,2T, ..., (M-l)T (2.9) 

We call the impulse response signal model, in the form of 
equation (2.9) , an implicit form of the new signal model. 

A simpler form of equation (2.9) can be either of equations 
(2.1j?) or (2.11) . 



N 



X(k) = E(k) + ^ A exp(S^k) 

n=l 


(2.10) 


N 


X(k) = 1 A^ (k) exp (S^k) 

n= l“ 


(2.11) 
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where A (k) = A , for k > T0 
n n — 

In equation (2.1), the residues are time-varying until 
the time at which the augmented function E(k) in equation 
(2.10) becomes zero. Thereafter the ~ ^n (^°^^stant) . 

To begin with equation (2.10), it is necessary to define 
a new independent variable as a parameter whose time behavior 
is not predictable along the positive time axis. It is con- 
venient to define this as a set of variables. By regarding 
them as independent, the equation (2.11) would be in the form 
of the equation (2.12) in the discrete digital data processing 
sense. 



assuming e (k) = E 6(k-n) 
n n 

The equation (2.14) is expressed in terms of a positive 
counting sequence in k. As in the equation, we might evaluate 
the equation (2.13) if we know the exact number of poles. But 
we know that it is impossible because we are in the region of 
the early time. So we write the equation (2.14) as a time 
shifted version such as the equation that has meaning from 
the observability viewpoint. 



N 

X(k) = E(k) + I A exp(S k) 
n=l ^ ^ 



( 2 . 12 ) 



T0-1 N 

X(k) = I e (k) + I \ exp (S^^k) 
n=l n=l 



(2.13) 



T0-1 N 

X(k) = I e^(k) + ^ A^ exp(S^k) 

n=l n=l 



(2.14) 
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where k=T0+m, T0+m-l, T0+l,T0,TJ2f-l, ... , Tj2,0 

for m ^ 0 

Here, the examples of the new signal model are presented 
in Appendix C, in the form of a decomposed signal. This syn- 
thetically generated signal is written in terms of 's 

and the sum of exponentials. In Chapter V, we present the 
decomposed signal from the synthetically generated data signal 
and that obtained by constructing the E(k) as the results 
of the non-linear parameter optimization processing. 

C. NOISE 

White Gaussian noise is assumed throughout this work. By 
introducing the noise into the new impulse signal model, we 
call the following equation the modeling function of a tran- 
sient response system and the complete and general form of 

the new signal model. 

oo N 

X(k) = N(k) + I e (k) + I A exp(S k) (2.15) 

n=l ^ n=l ^ ^ 
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III. PR0BLEI4S WITH THE TRADITIONAL METHODS 



In this chapter, a brief description of problems with 
the classic Prony's method, the Householder orthogonalization 
technique and the Eigenvalue method are provided. Each of 
these methods computes the location of poles in the s-plane 
but only the last two methods can estimate the number of poles 
prior to the discrimination of the poles. 



A. PRONY'S ALGORITHM 



Equation (2.3) is interpreted as the following difference 
equation assigning the values of the real part of the signal at 
specified sampling points in its left side of the expression. 



X{0) 



N 



I 

n=l 



A 



n 



X(k) 



N 



I 

n=l 



A 

n n 



(3.1) 



N 

X((M-l)k) = I A„Z^(M-l)k 

i, n n 
n=l 



And a polynomial equation (3.2) that has the same roots Z^ 

I a^z"" = 0 (3.2) 

n=l 
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c3.n bs combinsd with squation (3.1) to yisld ths Pirony's 
difference equation in the form of equation (3.3) 

N 

I a^X(n+k) = 0, for k=0,T,2T, ... , (M-1) T (3.3) 

n=l 

Then, we write equation (3.4) from equation (3.3) as 

N N-1 

I a X(n+k) = I a X(n+k) + a x(n+k) =0 (3.4) 

n=l “ n=l ^ ^ 

Using matrix notation, the equation (3.3) may be written as 
equation (3.5) or (3.6). 

N-1 

I a X(n+k) = -a^X(n+k) (3.5) 

n=l ^ 

Let a^^ = 1, then equation (3.5) becomes as 



^-1 ^-1 '^N 

where is N by N circulant matrix of sampled data 

is N by 1 Prony ' s coefficients matrix 
x^ is N by 1 row matrix of sampled data set 

^^'N+1' "* ' ^2N-1^ 



The Prony 's coefficient ' s are to be calculated by 
making use of the characteristics of the circulant matrix X. 

^-1 " ^^N-1 ^N-1^ ^-1 (3.7) 

As described in the above, this method can be applied after 

the number of poles are known to us, then the value of the 

S 's are to be found from the equation (3.8) . 
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(3.8) 



= LN(Z^)/k, for k > 2N 

B. HOUSEHOLDER ORTHOGONALIZATION METHOD 

A general expression of equation (3.2) may be represented 
by equation (3.9) assuming the number of poles are not known. 

^0 + ^1^^ + . . . + a^.Z^' = 0 (3.9) 

where N' is the unknown variable 
Then, the equation (3.9) is to be satisfied under the condi- 
tion that there are enough sampled data points m ^ 2N. 

where the N by 1 sampled data matrix is defined as 
x^ = [X(i), X(i+2) , . . . , x(M-N-i+l) ] 

Application of the successive orthogonalization through the 
Gram-Schmit process would produce the orthogonal vector set. 



0 - o^ , . . . , Ojj , } 

n-1 

where o = x - 7 <x ,o.> o. (3.11) 

n n . n' 1 i 
1=0 

If we set the to be 1, then the N by 1 x^ matrix is going 
to be in the form of equation (3.11) . 



n-1 



X = o + y<x,0. >o. 

n n n' 1 i 

1=0 



(3.12) 



So that the equation (3.12) holds as X=0M. In the above, 

M is the 2-D matrix of multiplication factors whose diagonal 
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elements are all I's. Here, an orthogonal vector o^ 
the corresponding data set to be orthogonal to all the 

previous vectors If we have had all the x vectors 

during i-1 times of successive orthogonalization process 
such that there is a non-zero orthogonal vector. We have to 
receive continuously the next data set check whether 

the orthogonal vector in the next step is zero or under a given 
threshold. If it vanishes, we can say that we have (i-1) 
poles. Our test against noise polluted data signal using the 
Householder method showed the inability of handling the early 
time signal as well as the data being heavily polluted. 

C. EIGENVALUE METHOD 

Again, from the equation (3.10), we can rewrite this 
equation as in the form of equation (3.14) . 

where X = • • • I’S)-! 

The eigenvector can be derived from the equation (3.15) by 
doing some mathematical manipulations of equation (3.14). 

^ N • ^ ^ ^ ^ N ' “N ' ^N • ' ^ 

From this equation, we find the eigenvalues correspond to 
eigenvectors to see if there is any eigenvalue approaching 
zero. If there is one zero or under threshold, we also 
consider the nximber of poles as N'-l. Even both the 
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Householder orthogonal! zation method and the Eigenvalue 
method can provide the means of calculating the number of 
poles, their algorithm (tending to fit the data points to 
the traditional signal model) showed the lack of generality. 
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IV. NON-LINEAR PARAMETER OPTIMIZATION APPROACH 



A. INTRODUCTION 

In 1963, Marquardt suggested an algorithm for the least- 
squares estimation of non-linear parameters [Ref. 6] when 
highly accurate parameter values are required. As we can 
see in the signal model in equation (2.14), that new signal 
model has multiple parameters that are functions of time. 
Here, we rewrite that equation again in the form of equation 
(4.1) . 

X(k) = N(k) + f e (k) + ^ ^n ®^P^S^k) (4.1) 

n=k " n=l 

where = 0t for k _>T0 

We are not interested in the value of itself, but 

include their estimations in order to contribute to the 
accuracy of S^^'s and A^'s. 

In other words, the extrapolation of the sum of exponen- 
tials, which is assumed to have all the poles, to points in 
the early time region may provide more accurate poles and 
residues simultaneously. 

The major advantage is the fact that we can make use of 
the high SNR data signal in the early time region with the 
basic information of poles that were derived in the late time 
region . 
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We will be processing on a data window whose vector 
length is to be increased one by one by moving the first 
element of that vector towa'rd the time-origin point. The 
optimization process evaluates the normal equation to find 
the local optimized set and finds the global optimized set 
as its final goal. 

Let us define an ERROR function as in equation (4.3) 

N-1 

ERROR{A^,S^,E(k) } = I [X’ (k) -X(k)]"^ (4.3) 

n=l 

where X* (k) is a measured data point 
To have minimized the least-squares error at every instant 
of measured time, the ERROR function has to be minimized by 
obtaining the global set of parameters being optimized. We 
also can write equation (4.3) in a more concrete manner as in 
equation (4.4). 

M-1 M-1 

ERROR {X} = I [X'(k) -X(k)]^= I e(k)^ (4.4) 

k=0 k=0 ^ 

where X = [E(k) ,A ,S ] 

Non-linear parameters such as must be obtimized in 

a way such that the more accurate values of have to be 

calculated as we increase the number of non-linear parameters 
e^. Users may define the numbers of parameter that 

are at least equal to or greater than the number of data 
points in the early time region. Now it is our task to find 
out the optimized value of and e^ through optimiza- 

tion processing, either in the global sense as it is represented 

in the equation (4.5) or using the normal equation (4.6). 
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9ERROR(X)/9X = jS 



(4.5) 



9e^(k)/9X(k) = 0 (4.6) 

where k=0,T,2T . . . , (M-l)K 

n=l,2, . . . , N 

With using the normal equation, as it is shown in the (4.6), 
a global set of non-linear set of parameters can be obtained 
through an iterative evaluation and fitting process. In our 
simulation work, the modified least-squares method was the 
basis of the equation (4.6). 

B. THE MODIFIED LEAST-SQUARES METHOD 

The error function is rewritten as in the expression of 
the equation (4.7). 

ERROR{k,x} = X’ (k) - X(k, (4.7) 

where X' (k) is the k-th sampled data 
X(k, x) is the model function 

is a vector containing the k-th sampled 

Let us define x^ to be an initial estimated value of x, then 

n 

a sequence of approximations to the optimized value is to be 
generated by the equation (4.8). 

m+1 ^ ^ ^ jT j-1 jT ERROR (x"') (4.8) 

*•111111 mm m 

where J is the numerical Jacobian matrix evaluated 
m 

m is iteration number of successive optimization 
D is a diagonal matrix equal to the diagonal of 
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T 

J J 
m m 



is a Marquardt parameter 

The number of total iterations can be controlled by the 
threshold which is defined as in the equation (4.9). 

. m^* 1 m ! A Q \ 

d = X - X (4.9) 
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V. TEST RESULTS AND PERFORMANCE EVALUATION 



A. INTRODUCTION 

In order to establish the ability of the program listed 
in Appendix B to extract the poles and residues correctly, 
three different simulated signals, each of which is polluted 
with infinite, 30 dB and 15 dB SNR, were created by the syn- 
thetic signal data generation routine. These 3 sets of signal 
data were chosen to span a range of possible situations. 

In the context of the transient signals and additive sta- 
tionary Gaussian noise being used , the definition of SNR is 
in terms of a ratio of energy quantities intergrated over the 
entire 20 nsec time window. 

TABLE I 

SIMULATED SIGNAL 1 
POLES 

RESIDUES ( Nep. GHZ ) 

1.0 + il.0 -1.0 + il.0 

1.0 - 11.0 -1.0 - 11.0 

Simulated signal 2 consists of 2 sets of pairs of complex 

conjugate poles and residues, which is an extrapolation of the 
simulated signal 1, using its parameters. Simulated signal 3 
has 3 sets of pole-residue pairs, extrapolating from signal 2. 

There are 5 options in the signal generation program in 

choosing a particular function of E(k) . In this simulation 

work, a trapezoidal wave form of E(k) was used with T0=0 .42 

So that the unknown parameters that are to be augmented 
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nsec . 



at every processing step will be a maximum of 10 plus 4 times 
the number of poles at the final processing stage, when we 
have 512 data points within a 20 nsec time window. 

TABLE II 

SIMULATED SIGNAL 2 



RESIDUES 


POLES 

(Nep. GHZ) 


1.0 + il.0 


-1.0 + il.0 


i.0 - il.0 


-1.0 - il.0 


0.5 + i0.5 


-2.0 + i2.0 


0.5 + i0.5 


-2.0 - j2.0 



TABLE III 



SIMULATED 


SIGNAL 3 


RESIDUES 


POLES 

(Nep. GHZ) 


1.0 + il.0 


-1.0 + il.0 


1.0 - il.0 


-1.0 - il.0 


0.5 + i0.5 


-2.0 + i2.0 


0.5 - i0.5 


-2.0 - 12.0 


0.25 + 10.25 


-3.0 + i3.0 


0.25 - 10.25 


-3.0 - i3.0 


The results of varying the 


number of additional data 



points from zero (all the data points are from late time 
region) to ten extra data points (10 extra points are from 
early time region are added to the points of late time 
region) are contained in Table IV through XII. These nine 
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cases correspond to the 3 synthetic signals, each having 
3 additive noise levels. The 3 synthetic signals are plotted 
in Appendix C. In the next section, the accuracy of pole 
extractions, as indicated in the tables are displayed by way 
of graphical pole maps. In addition, reconstructed waveform 
obtained from the parameter extractions are compared graphi- 
cally to the original waveforms. 
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TABLE IV 



PARAMETER OPTIMIZATION FOR SIGNAL 1 (NOISE FREE) 



B, 



RESULTS 



TARGET TYPE:TGT-1 
**AVEfCRM TYFE:hCNFTR 
CONTACT LATEiCEC 15 
FILS MJH£:FILE002 

NUMB. CF pole: 2 

TABLE CF RESIDUES AND POLES 



PAIR # 
1 
2 



RES-REAL 
1. OCCGOOOO 
1. CCCCOOCO 



RES-IMAG POLE-REAL 
1.00000000 -l.OOOOCOOO 
-1. CCOCCOOO -l.COOCCOOG 



FOLE-IMAG 

1.00000000 

- 1.00000000 



RESULTS CF OFIIMZATICN hITh 0 EXTRA DATA PCINTS 



PAIR 


U- 


RES. -REAL 


RES.-IMAG 


POLE-REAL 


FOLE-IMAG 


1 




C. SSS96610 


1. 00001592 


-0.99994705 


C. 99998971 


2 




0.95S96563 


-1. 00001375 


-0.99994777 


-0.99998971 


RESULTS 


CF GPTIMIZATICN WITH 2 


EXTRA DATA PCINTS 


PAIR 


a 


RES. -REAL 


RES.-IMAG 


FCLE-REAL 


FOLE-IMAG 


1 




1. 0CC0O251 


0. 99999667 


-1 .G00CC025 


1.0C000085 


2 




1. CCC00251 


-0. 99999667 


-1 .00000026 


-1.0C000085 


RESULTS 


CF OFTIMIZATICN WITH 4 


EXTRA DATA PCINTS 


PAIR 




RES. -REAL 


RES.-IMAG 


POLE-REAL 


FOLE-IMAG 


i 




1.0CC013C8 


1. CC000329 


-1 .C00C1536 


0.99999999 


2 




1. 0CCC1311 


-1. COOOC330 


-1 .00001545 


-1.00000002 


RESULTS 


CF OFIIMZATICN WITH 6 


EXTRA DATA PCINTS 


PAIR 




RES. -REAL 


RES.-IMAG 


POLE-REAL 


FOLE-IMAG 


1 




C.9SS25213 


1. 00003040 


-0.99979429 


0.99998208 


2 




0. 9SS85210 


-1. CC003047 


-0.99979441 


-0.99993213 


RESULTS 


CF OFTIMZATICN WITH 8 


EXTRA DATA PCINTS 


PAIR 


4 


RES .-REAL 


RES.-IMAG 


POLE-REAL 


FOLE-IMAG 


1 




1. CCCCC036 


1. GCOOlOll 


-1.00000472 


0.95999818 


2 




1. OCCOC026 


-1.0C001041 


-1.C00CC54C 


-0.99999838 


RESULTS 


GF OFTIMIZATICN WITH 10 


EXTRA DATA PCINTS 


PAIR 


n 


RES. -REAL 


RES.-IMAG 


POLE-REAL 


FOLE-IMAG 


1 




1.0CC00457 


C. 99959062 


-0.99999674 


1.C0000165 


2 




1 • C C C 0 0^ ^3 


-0. 99999067 


-0 .99999735 


-1 .00000168 
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TABLE V 



PARAMETER OPTIMIZATION FOR SIGNAL 1 
(SNR=30 d3) 



TARGET TYP£;TGT-1 
WAVEFORM TYPE :l-CN30TR 
CONTACT CATE:CEC 15 
FILE NUME;fILE003 

NUMB. CF POLE; 2 





TABLE 


CF RESIDUES 


AND POLES 






PAIR 


# RES-REAL 


RES-IMAG 


POLE-REAL 


POL 5- I M 


ACj 


1 


l.OCCCOOCO 


1. 00000000 


-l.OOOCQOOO 


1 .OOOOOo 


00 


2 


1. GCCCOOCO 


-1. CCOCCCOO 


-i.COCCCOOC 


-l.OOQCOO 


09 


RESU 


LTS CF OPTIMIZATICN WITH 0 


EXTRA DATA PCINTS 




PAIR 


f# RES. -REAL 


FES.-IMAG POLE-REAL 


POLE-IM 


AG 


1 


1.01CC6572 


0.92346346 -0.77676493 


1.003700 


76 


2 


1.03315188 


-1.10738236 -1.30247523 


0.991974 


36 


RESULTS CF 0P1IMIZATICN i^ITH 2 


EXTRA DATA PCINTS 




PAIR 


# RES. -PEAL 


RES.-IMAG 


POLE-REAL 


POLE-IM 


AG 


1 


1 .0107C336 


1. CC36C284 


-1.00951167 


0.999793 


04 


2 


1 .0107C339 


-1. 00360275 


-1 .00951163 


0.999793 


06 


RESULTS CF OPTIMIZATICN WITH 4 


EXTRA DATA PCINTS 




PAIR 


ft RES. -REAL 


RES.-IMAG 


POLE-REAL 


POLE-IM 


AG 


1 


1 .01066217 


1.003o7655 


-1.00953849 


0.999777 


51 


2 


1.01066632 


-1.0C367409 


-1 .00953782 


-0.999777 


35 


RESULTS CF OPTIMIZATICN WITH 6 


EXTRA DATA PCINTS 




PAIR 


ft RES. -REAL • 


RES.-IMAG 


POLE-REAL 


POLE-IM 


AG 


1 


1 .01059740 


1.00357123 


-1 .00930702 


0.999785 


80 


2 


1.01059679 


-1.00357126 


-1.C0930751 


0.999735 


74 


RESULTS CF OPTIMIZATICN WITH 8 


EXTRA DATA PCINTS 




PAIR 


# RES. -PEAL 


RES .-IMAG 


POLE-REAL 


POLE-IM 


AG 


I 


1. 010679^0 


1.00362571 


-1.00950587 


0.999786 


86 


2 


l.C106£C33 


1.00362636 


-1.C095058E 


-0.999786 


99 


RESULTS CF OPTIMIZATICN WITH 10 


EXTRA DATA PCINTS 




PAIR 


ft RES. -PEAL 


RES .-IMAG 


PCLE-REAL 


POLE-IM 


AG 


1 


1 .01069621 


1.00360379 


-1 .C095C576 


0.999791 


75 


2 


1.01C69o09 - 


1 .00360399 


-1 .00950630 


-0 .999791 


77 



31 



TABLE VI 



PARAMETER OPTIMIZATION FOR SIGNAL 1 
(SNR=15 dB) 



TARGET 

WAVEFORM 

CONTACT 

FILE 

NUM6. OF 


TYPE:TGT-1 
TYPE :hDN15TR 
CATE:CEC 20 
NUME :FILE004 
POLE: 2 










TABLE 


CF RESIDUES 


AND POLES 




PAIR » 
1 
2 


REB-REAL 
1. CCCCCOCO 
l.CCCOOOOO 


RES-IMAG 
1. CCOOOOOO 
-1. COOOOOOO 


POLE-REAL 
-1.00000000 
-l.OOOCCOOC - 


POLE-IMAG 
1 .00000000 
1.00000000 


RESULTS 


CF OPTIMIZATICN WITH 0 


EXTRA DATA PCINTS 


PAIR # 
1 
2 


RES. -REAL 
0. 56CA79E2 
1.59727710 


RES.-IMAG 
1. 29922892 
-0. 80336947 


PCLE-REAL 
-0.5789C084 
-1.57012329 - 


POLE-IMAG 
0.968863 04 
0.914193 74 


RESULTS 


CF OPTIMIZATICN WITH 2 


EXTRA DATA PCINTS 


PAIR # 
1 
2 


RES. -REAL 
1.06C57099 
1. 0EC57923 


RES.-IMAG 
1. 02070489 
-1. 02071795 


POLE-REAL 
-1.05306349 
-1.05306276 - 


POLE-IMAG 

0.99E69879 

•0.99869878 


RESULTS 


CF OPTIMIZATICN WITH 4 


EXTRA DATA PCINTS 


PAIR # 
1 
2 


RES. -REAL 
1.06C53386 
1. 06C52818 


RES.-IMAG 
1. 02072970 
-1. 02073354 


POLE-REAL 
-1.05304404 
-1.052C4091 - 


POLE-IMAG 

0.99868990 

0.99869088 


RESULTS 


CF OPTIMIZATICN WITH 6 


EXTRA DATA PCINTS 


PAIR # 
1 
2 


RES .-REAL 
1. 0EC51499 
1. 08054674 


RES.-IMAG 
1. 020 73547 
-1. 02072794 


POLE-REAL 
-1.05295266 
-1.05313212 - 


POLE-IMAG 

0.99873653 

0.99864419 


RESULTS 


CF OPTIMIZATICN WITH 8 


EXTRA DATA PCINTS 


PAIR 4 
1 
2 


RES. -REAL 
1. 02A49720 
1. 07C48258 


RES.-IMAG 
1. 09889916 
-1. 10030318 


pole-real 

-0.66485231 
-1.63907420 - 


POLE-IMAG 

1.01775339 

0.93004124 


RESULTS 


CF OPTIMIZATION WITH 10 


EXTRA DATA PCINTS 


PAIR n 
1 
2 


RES. -REAL 
1.02166928 
1. 02251816 


RES.-IMAG 
1. 11984550 
-1. 07048377 


POcE-REAL 
-1. 566753 Ir 
-0.66239621 - 


POLE-IMAG 
C. 92761091 
1.02160947 
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TABLE VII 



FOLt-IM AG 
1.00000000 
- 1.00000000 
2.00000000 
- 2.00000000 



RESULTS CF CFliMZATICN 0 EXTRA DATA PCJNTS 



PAIR 




RES .-BEAL 


RES.-IMAG 


POLE-REAL 


FOLE-IMAG 


1 




1. 0CCCCC54 


0. 99998637 


-1.00000054 


1.00000162 


2 




1.GCCC0019 


-0.99998792 


-0.99999323 


-1.000001 80 


3 




C. 49S93655 


0. 50005618 


-2.0000551 8 


1.99995365 


4 




0. 499S5131 


-C. 50004431 


-2.00024952 


-1.999941 15 


RESULTS 


CF OFTIRIZATICN WITH 2 


EXTRA DATA POINTS 


PAIR 




RES .-REAL 


RES.-IMAG 


PCLE-REAL 


FOLE-IMAG 


1 




C. 5SSC6513 


C. 99549512 


-0.99792594 


0.999981 85 


2 




0.9S9064C5 


-C. 99549467 


-0.99792606 


-0.999987 94 


3 




0.47735342 


0.48850313 


-1.95118748 


1.99896225 


4 




0. 47735838 


-C. 48850442 


-1.953C6779 


-1.99350918 


RESULTS 


CF OFTIPIZATICN WITH 4 


EXTRA DATA POINTS 


PAIR 




RES. -REAL 


RES.-IMAG 


PCLE-REAL 


FOLE-IMAG 


1 




C. 999Cc256 


0. 99549224 


-0.99793418 


0.99998541 


2 




0.99906321 


-C. 99549228 


-0.99791438 


-0.99998413 


3 




0.47737187 


0. 48348088 


-1.95232122 


1.99855642 


4 




0.47 736975 


-C. 48847526 


-1.95190026 


-1.99893520 


RESULTS 


CF 0P7IMIZATICN WITH 6 


EXTRA DATA POINTS 


PAIR 


# 


RES. -REAL 


RES.-IMAG 


PCLE-REAL 


FOLE-IMAG 


1 




C. 999C6743 


C. 99549467 


-0.99795003 


0.99998501 


2 




0.999C6717 


-0. 99549457 


-0.99792354 


-C.99993468 


3 




C. 47739612 


0. 48847662 


-1.95179495 


1.99874576 


4 




0.47739443 


-C. 48847686 


-1.95248987 


-1.99875751 


RESULTS 


CF CFTIMZATICN WITH 8 


EXTRA DATA POINTS 


PAIR 




RES. -REAL 


RES.-IMAG 


POLE-REAL 


FOL E— IM AG 


1 




i.01£14eC9 


C. 98882407 


-1.12954971 


1.01650809 


2 




1.0C949031 


-0.98875294 


-0.88470066 


-0.99067696 


3 




0.7E175371 


0. 36507609 


-2.03189450 


2.04 7446 71 


4 




C. 14319333 


-0. 56814961 


-1.95892411 


-1.92582259 


RESULTS 


CF GPTIPIZATICN WITH 10 


EXTRA DATA PCI^ 


ITS 


PAIR 




RES. -REAL 


RES.-IMAG 


POLE-REAL 


FOLE-IMAG 


1 




1. 02706396 


0. 98982102 


-1.16328021 


1.01920699 


2 




1.02C53154 


-C. 99088746 


-0.8681587C 


-0.99077552 


3 




0.34656670 


0. 59967646 


-1.95544269 


1.95297875 


4 




0. 615 55028 


-0. 37246541 


-2.00784742 


-2.04312205 



PARAMETER OPTIMIZATION FOR SIGNAL 2 
(NOISE FREE) 

TARGET TYPE:TGT-2 
WAVEFORM TYPE:hCNFTR 
CONTACT CATE:CEC 15 
FILE AUM£:fILE002 

NJMB. CF POLE: 4 

TAELE CF RESIDUES AND POLES 



PAIR # RES-REAL 

1 l.OCCCOCCO 

2 l.GCCOOOOO 

3 0.5CCOCOCO 

A C.5CCGOQCO 



RES-IMAG 
1. COOOOOOO 
•l.CCOOOOOO 
0. 5C00C000 
•0. 5C0CC000 



POLE-REAL 

- 1.00000000 

-l.OOOCCOOO 

- 2.00000000 

-2.000CC00C 
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TABLE VIII 



PARAMETER OPTIMIZATION FOR SIGNAL 2 
(SNR=30 dB) 

TARGET 7YPE:7GT-2 
WAVEFORM TYPE :h0N30TR 
CONTACT DATE:CEC 15 
FILE NUMB:FILEQ03 

NUMB, CF POLE: 4 

TABLE CF RESIDUES AND POLES 



PAIR « RES-REAL 

1 l.OCCQOOOO 

2 1,CCCCC0C0 

3 C. 5CCOCOOO 

4 C.5CCCC0C0 



RES-IMAG 
1. COOCOOOO 
•1. CCOOOOOO 
0. 5COCCDOO 
•0. 50000000 



pole-real 

-I .oocccooo 

-i .00000000 
-2.000CC00C 
- 2 .00000000 



FCLE-IMAG 
1.00000000 
- 1 .00000000 
2 .00 C 000 CO 
-2 .00003000 



RESULTS CF 0FTIMI2ATICN WITH 0 EXTRA DATA POINTS 



PAIR 




RES. -REAL 


RES.-IMAG 


FCLE-REAL 




1 




C. 99412945 


0. 99056805 


-0.92918755 


0 


2 




C. 99*48260 


-C. 99069937 


-i .06493363 


-1 


3 




C. CCCGOOlO 


C. 56023544 


-1.96891695 


1, 


4 




0. 87406384 


-0. 32304349 


-2.04566871 


-2 


RESULTS 


CF OFTIMIZATICN WITH 2 EXTRA DATA PCINTS 


PAIR 




RES. -REAL 


RES.-IMAG 


FOLE-REAL 


1 


1 




0.99906513 


0. 99549512 


-0.99792594 


o' 


2 




0. 9 99C64C5 


-0. 99549467 


-0.99792606 


-0. 


3 




0. 47735342 


0. 48850313 


-1.95118748 


1. 


4 




0. 47735838 


-0. 48850442 


-1.95306779 


- i ^ 


RESULTS 


CF OFTIMIZATICN WITH 4 EXTRA DATA PCINTS 


PAIR 




RES. -REAL 


RES.-IMAG FOLE-REAL 


1 


1 




0.99906256 


0. 99549224 


-0.99793416 


0 


2 




C. 99906321 


-0. 99549228 


-0.99791438 


-0. 


3 




0.47737187 


0. 48848088 


-1.95232122 


1. 


4 




0.47736977 


-0. 48847526 


-1.95190026 


-1, 


RESULTS 


CF OFTIMIZATICN WITH 6 EXTRA DATA PCINTS 


PAIR 


# 


RES. -REAL 


RES.-IMAG 


FOLE-REAL 




1 




0.99906743 


C. 99549467 


-0.99793003 


0 


2 




0. 99906717 


-C. 99549457 


-0.99792354 


-0 


3 




0. 47739612 


0. 48847662 


-1.95179495 


1. 


4 




0.47739443 


-0. 48847686 


-1 .9524898? 


-1 


RESULTS 


CF OFTIMIZATICN WITH 8 EXTRA DATA POINTS 


PAIR 




RES. -real 


RES.-IMAG 


FOLE-REAL 




1 




1.015148C9 


C. 98882407 


-1 .1295497 1 


1 


2 




1.CC949031 


-C. 98875294 


-0.88470066 


-0 


3 




0.76175371 


0. 36507609 


-2.0318945C 


2 


4 




0. 14319333 


-0. 56814961 


-1 .9589241 1 


— 1 


RESULTS 


CF OFTIMIZATICN WITH 10 


EXTRA DATA PCINTS 


PAIR 




RES .-REAL 


RES.-IMAG 


FOLE-REAL 


1 




1.027C6396 


0. 98982102 


-1 .16328021 


1 


2 




1. 02C53154 


-C. 99088746 


-0.86815870 


-c 


3 




0. 34656670 


0. 59967646 


-1 . 95544269 


1 


4 




0.61555028 


-0. 37246541 


-2 .00784742 


-2 



FOLE-IMAG 



FOLE-IMAG 



POLE-IMAG 

.SS998541 

.99998413 

.99855042 



PCLE-IMAO 



POLE— I M AG 



FOLE-IMAG 
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TABLE IX 



PARAMETER OPTIMIZATION FOR SIGNAL 2 
(SNR=15 dB) 

TARGET TYPE:7GT-2 
WAVEFGRM TYPE:hON15TR 
CONTACT CATE:CEC 15 
FILE AUM£:F1LE004 

NUMB, CF POLE: 4 

TABLE CF RESIDUES AND POLES 



PAIR )f 
1 
2 

3 

4 



RES-REAL 
l.OCCCOOOO 
1, CCCJQOOO 
C. 5CCOCOOO 
C. 5CCCCCC0 



RES-IMAG 
1 . 00000000 
•i. CCOOOOOO 
0. 50000000 
•0. 50000000 



POLE-REAL 

- 1,00000000 

-1,OOCCOOOO 

-2.000CC0OC 

- 2,00000000 



POLE-IMAG 

1,00000000 

- 1,00000000 

2,00000000 

- 2,00000000 



RESULTS CF OFTIRIZATICN UlTh 0 EXTRA DATA POINTS 



PAIR 




RES .-REAL 


RES.-IMAG 


POLE-REAL 




1 




C.SBC91205 


0. 96731225 


-0.97658042 


0 


2 




C. 9 6CS26C3 


-0. 56731915 


-0.57544458 


-0 


5 




0.52307777 


0. 22152280 


-1.54313287 


2 


4 




C. CCCOCOlO 


-0.42428522 


-1 ,38794994 


-1 


RESULTS 


CF OFTIMIZATICN l<ITh 2 


EXTRA DATA POINTS 


PAIR 




RES. -REAL 


RES.-IMAG 


POLE-REAL 




1 




C.S£51i311 


0. 96627888 


-0.59350301 


1 


2 




0. 9A526467 


-0. 96652064 


-0.95966500 


-0 


3 




O.GCCOOOlO 


0. 49325107 


-1 .46SC5167 


1 


4 




0. 5545C498 


-0. 13648901 


-1.66572130 


-2 


RESULTS 


CF OFIIMIZATICN WITH 4 


EXTRA DATA POINTS 


PAIR 




RES , -REAL 


RES.-IMAG 


POLE-REAL 




1 




C. 56128516 


0. 56675315 


-0.98935765 


0 


2 




0. 56144733 


-0. 56682559 


-0.56226170 


-0 


3 




O.OCCCOOlO 


0. 54101555 


-1 .50065226 


1 


4 




0. 5C2513B3 


-0. 09527370 


-1 .60980060 


-2 


RESULTS 


CF OPTIHIZATICN ViITh 6 


EXTRA DATA POINTS 


PAIR 


it 


RES. -REAL 


RES.-IMAG 


POLE-REAL 




1 




1.06254479 


0. 93763540 


-1.25945424 


1 


2 




1. C 1484084 


-0. 54055523 


-0.7727797C 


-0 


3 




0.62674427 


0. 16116490 


-1 .75036751 


2 


4 




C. CCCOOOlO 


-0. 5197C056 


-1 .59036824 


-1 


RESULTS 


CF OPTIMIZATICN k<ITh 8 


EXTRA DATA POINTS 


PAIR 




RES. -REAL 


RES.-IMAG 


POLE-REAL 




1 




0.563C4511 


0. 56782052 


-0.58445408 


0 


2 




0, 56312358 


-0. 56782270 


-0.5654CS93 


-0 


3 




O.OCCOOOlO 


0.35715508 


-1.36627593 


1 


4 




0. 54632361 


-0. 25136131 


-1 .56644065 


-2 


RESULTS 


CF OPTINIZAT ICN l»ITh 10 


EXTRA DATA PCINTS 


PAIR 


it 


RES. -REAL 


RES.-IMAG 


POLE-REAL 




1 




1. 0C571067 


0. 56985420 


-1 .20987398 


1 


2 




1. 0C285344 


-C. 57024550 


-0.81824047 


-0 


3 




0, 576C7620 


0. 19510021 


-1.66621598 


2 


4 




C. OCCCOCIO 


-0. 50028280 


-1.53334141 


-1 



PGLE-IMAG 



POLE-IMAG 



POLE-IMAG 



POLE-IMAG 



POLE-IMAG 



POLE-IMAG 
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TABLE X 



PARAMETER OPTIMIZATION FOR SIGNAL 3 
(NOISE FREE) 



TARGET 
WAVS FCBM 

Contact 

FILE 

NUM8. CF 



TYPE :7GT-3 
TYPE :hCNFTR 
CATE :C£C 15 
NUMB :FIL£00A 
POLE: fc 



TABLE CF RESIDUES AND POlES 



PAIR # RES-REAL 

1 1.0CC30a03 

2 l.CCCCCOCO 

5 0.5CCQC000 

A C.5CCCC0C0 

5 C.25CGG0C0 

6 C.25CCCGC0 



1 . 

0 . 

- 0 . 

0 . 



RES- 1 MAG 
OOOOCOOO 

ccocccoo 

5CD00000 

5C0CC000 

2500G000 



-0. 25000000 



FCLE-REAL 
-1 .OOOOCOOO 
-i.OOCCCOOC 
-2.000CC00C 
- 2.00000000 
-3. COOCCCOO 
-3 .00000000 



FOLE-IMAG 
1 .000000 00 
- 1.00000000 
2.00000000 
- 2.00000000 
3.00000000 
-3.00000000 



RESULTS CF CrTiMiZATICN UlTH 0 EXTRA DATA POINTS 



PAIR # RES. -REAL 

1 1.0CC00029 

2 0.9CC99965 

3 C.5CC0C021 

4 • 0.49999957 

5 C.25CCCC49 

6 0.24999763 



RES. -I MAG 
0. 99999951 
-1. 00000047 
0. 49999977 
-0. 50000051 
C. 25000087 
-0. 24999958 



FCLt-REAL 
-1 .C00CC012 
-0.99999985 
-1.99999733 
-2 .00000290 
-2.99676241 
-3.00125324 



FOLE-IMAG 
1 .00000000 
-0.99999999 
1.99999965 
-2.00000029 
2.99999506 
-2.99999619 



RESULTS CF CFTIMIZATICN WITH 2 EXTRA DATA PCINTS 



PAIR H RES. -REAL 

1 1.0CCQ0033 

2 0.99999970 

3 C.49999897 

4 C.5CC0C019 

5 0. 24999312 

6 0. 25C00768 



RES.-IMAG 
1. CCCOC006 
•0. 99999990 
C. 49999976 
•0. 49999970 
C. 25000063 
'C. 24999482 



FCLE-REAL 
-0.99999640 
-1 .00000359 
-2.00000082 
-1.99999827 
-2.99929864 
-3.0006786B 



FOLE-IMAG 

0.99999972 

-1.00000028 

.1.99999976 

-2.CC000017 

2.99995040 

-5.00005033 



RESULTS CF OFTIMIZATICN WITH 4 EXTRA DATA PCINTS 



PAIR # 
1 
2 

3 

4 

5 

6 



RES. -REAL 
1. 0CC00031 
1. CCCC00C8 
C.5CC0CC29 
0. 5CC0C073 
0. 24999713 
C. 25C003C6 



RES.-IMAG 
1. CC000004 
-0. 99999999 
0. 4S999878 
-0. 49999913 
0. 250 00137 
-0. 24999453 



FCLE-REAL 
-0. 99999929 
-1 .00000084 
-2.C0000158 
-1.9999974E 
-3.00006822 
-2.99992441 



FOLE-IMAG 
0.99999999 
- 1.00000002 
2.00000057 
-1.999999 86 
3.00000834 
-2.99999365 
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RESULTS CF CFTIMZATICN WITH 6 EXTRA DATA PCINTS 



PAIR # RES. -REAL 

1 C. SSS99967 

2 0.9^999949 

3 C.5CCOC350 

A C.5CCQC357 

5 0.25C02351 

6 0.24999C17 



RES.-IMAG 
1. CC000006 
>1. 00000005 
0. 49999519 
•C. 49999515 
C. 25003725 
•C. 25002880 



PCLE-REAL 
“0.99999978 
-1.00GCC003 
-2.000CC41 1 
-2.C00CC5o5 
-3.01254931 
-2.9S71C6S7 



FOLE-iM AG 
0.99999979 
-1.0CC00016 
2.000001 16 
-2.000001 16 
2 .99 9940 2c 
-2.99957756 



RESULTS CF OFTIMZATICN UITh 8 EXTRA DATA PCINTS 



PAIR # RES. -REAL 

1 i.occaoooo 

2 i.OCCQOOOl 

3 0.5CCCC063 

4 0. 49999755 

5 0.25C667S7 

6 0. 24913768 



RES.-IMAG 
0. 99999980 
-1. COOCCOlO 
0.49995686 
-C. 50003918 
C. 25234C01 
-0. 24765392 



pole-real 

-0.99999966 
-1 .000CC03 1 

-2 .oooccesi 

-1 .99998120 
-3.00656519 
-2 .99334543 



FOLE-IMAG 
1. 000000 2o 
-0.999999 75 
1.99969134 
-2.00030891 
3.00052843 
-2.99947896 



RESULTS CF OPTINIZATICN UITh 10 EXTRA DATA PCINTS 



PAIR # RES. -REAL 

1 1.QCC00004 

2 0.99999994 

3 C.5CC00006 

4 C.5CC0CC19 

5 0.25CC0C13 

6 C.25C0CC24 



RES.-IMAG 
0. 99999995 
-1. CC000008 
0. 50000023 
-0. 5C00C011 
0. 25000022 
-C. 25000077 



POLE-REAL 
-1 .00000147 
-0.99999854 
-2.000C1924 
-1.99998131 
-2.99992168 
-3.0000816 1 



FOLE-IMAG 
1.00000006 
-0.999999 94 
1 .599998 46 
-2.00000154 
2.999975 83 
-3.000023 82 
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TABLE XI 



PARAt'lETER OPTIMIZATION FOR SIGNAL 3 
(SNR=30 dB) 



TARGET 
WAVEFORM 
CONTACT 
F IL£ 

NUMB. OF 



TYP E:TGT-3 
TYP5;HON30TR 
OATE.'OEC 15 
NJMB: FI LE004 
POLE: 6 



TABLE OF RESIDUES AND POLES 



PAIR 

1 

2 

3 

4 

5 

6 



1 . 
1 . 
0 . 
0 . 
0 . 
0 . 



RES-REAL 

00000000 

00000000 

50000000 

50000000 

25OOOC0O 

25000000 



RES- I MAS 
I .00000030 
- 1 .00000000 
0 .50000000 
-0.50000000 
0 .25000000 
-0 .25000000 



POLE-REAL 
- 1 .00000000 
- 1.00000000 
- 2.00000000 
-2 .00000000 
- 3. 00000000 
-3.00000000 



POLE-IMAG 
1.00000000 
- 1.00000000 
2. 00000000 
- 2.00000000 
3. 00000000 
-3.00000000 



RESULTS OF OPTIMIZATION WITH 0 EXTRA DATA POINTS 



PAIR U RES. -REAL 
1 0.99145389 

■ 2 0. 99145320 

3 0.42807658 

4 0.42846454 

5 0.11656813 

6 0.12939770 



RES. -I MAS 
0 .98947089 
-0 .98947086 
0. 52733218 
-0 .52679082 
0 .23 647261 
-0 .21591933 



POLE-REAL 
-0.99318094 
-0 .99310553 
-1.99204259 
-1 .92072086 
-2.27124794 
-2.24119798 



POLE-IMAG 
1.00001443 
-1.00001237 
1. 97341708 
- 1.99797545 
2.83611689 
-3. 06674116 



RESULTS OF OPTIMIZATION WITH 2 EXTRA DATA POINTS 



PAIR 

1 

2 

3 

4 

5 

6 



RES. -REAL 
0. 99147340 
0.99147541 
0.42927584 
0.4 29 3 8 7 21 
0 .13113198 
0. 10954599 



RES .- IMAG 
0.98938441 
-0 .98938831 
0 .52783363 
-0 . 52733004 
0 .20433430 
-0.24676563 



POLE-REAL 
-0.99323692 
-0.99299518 
-1.96123606 
-1.96042724 
-2.22870752 
-2 .26661642 



POLE-IMAG 
1.000021 18 
-1.00001483 
1.86208680 
-1.98830969 
3.07350223 
-2.8335683 1 



RESULTS OF OPTIMIZATION WITH 4 EXTRA DATA POINTS 



PAIR H. 
1 
2 

3 

4 

5 

6 



RES. -REAL 
0. 99149793 
0. 99149653 
0.42898254 
0.42897978 
0. 1 3728 153 
0 .1 0518882 



RES .- IMAS 
0 .98938042 
-0. 98937734 
0 .52744729 
-0. 52746013 



POLE-REAL 
-0.99300311 
-0.99323733 
-1.95708423 
-1 .96340652 



0.20601283 -2.24720729 

-0.24133713 -2.27543055 -2.83274078 



POLE-IMAG 
1.00001643 
-1.000 02228 
1.98692609 
-1.98537735 
3.07453032 
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RESJLTS OF OPTIMIZATION WITH 6 EXTRA DATA POINTS 



PAIR 4 RES. -REAL 

1 D.99L8L50& 

2 0.9918L556 

3 0.4 30 71914 

4 0.43071211 

5 0.15055C56 

6 0.10533241 



RES. -I MAO 
0 .98944282 
-0 . 98 9‘t4359 
0 . 5243 06o5 
-0 .52479742 
0.25270679 
-0 .20048107 



POLE-REAL 
-0.99314260 
-0.99332478 
- 1 .97634267 
-1 .93934530 
-2.28536097 
-2.19638655 



POLE-IMAG 
1. 00003060 

- 1.00003549 
1. 98151201 

- 1.99173878 
2.85503cl9 

-3. 07720611 



RESULTS OF OPTIMIZATION WITH 8 EXTRA DATA POINTS 



PAIR 

1 

2 

3 

4 

5 

6 



RES. -REAL 
o. 98171 120 
0 .981 71 167 
0.3 0805752 
0.21532625 
0 . 11572596 
0. 0 7739011 



RES .-IMA3 
0.93981746 
-0 .98981744 
0 .H9340155 
-0. 51226642 
0 .13443903 
-0.19217344 



POLE-REAL 
-0.99074872 
-0. 99060276 
-1 .38871252 
-1.56132982 
-1 .99741621 
-1.89148931 



POLE -I MAG 
0.999472 32 
-0. 99949877 
1.87285997 
-2.021278 1 1 
3. 12706630 
-2.85558132 



RESULTS OF OPTIMIZATION WITH 10 EXTRA DATA POINTS 



PAIR 4 RES. -REAL 

1 0.99337417 

2 0.99337 150 

3 0.442 3 9 942 

4 0. 4423 9753 

5 0.15478007 

6 0.15377770 



RES .-IMAG 
0.99015792 
-0.99015803 
0 .52004030 
-0. 52004763 
0 .23652003 
-0.23305063 



POLE-REAL 
-0.99418591 
-0.99369290 
-1 .95770753 
-1 .96547995 
-2.41823246 
-2 .48497797 



POLE-IMAG 
1.00008890 
-1.00007652 
1.991 13003 
-1.9874627 1 
2. 8709955 1 
-3.065 15871 
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TABLE XII 



PARAMETER OPTIMIZATION FOR SIGNAL 3 
(SNR=15 dB) 



TARGET 
WAVE FORM 
CONTACT 
FILE 

NUMB. CF 
PAIR U 


TYPE ;TGT-3 
TYPE :FCM5TR 
DATE :CEC 15 
NUME : FILE004 
POLE: t 

TA^E_CF_RESIDUES 
RES-REAL ficS-IMAG 


AND POLES 
POLE-REAL 


FOLE-Ii-i 


Au 


1 


l.OCCCOOOO 


1. 00000000 


-l.OOOCCOOO 


1 .000000 


00 


2 


l.CCCOOOOO 


-1. ccoccooo 


-i.OOOCCCOC 


-1 .000030 


00 


3 


C.5CCOCQOO 


0. 50000000 


-2 .OOOCCOOQ 


2.000000 


00 


4 


C.5CCCC0C0 


-C. 5COOOOOO 


-2.00000000 


-2.000000 


00 


5 


C.EfCCOOCO 


C. 25000000 


-3.000CCOOC 


3.000000 


go 


6 


C.25CGC000 


-0. 25000000 


-3 .00000000 


-3.000000 


00 



RESULTS 


CF OPTIMZAT 


ICN WITH 0 


EXTRA DATA POINTS 


PAIR 


# 


RES .-REAL 


RES.-IMAG 


POLE-REAL 


FGLE-IMA& 


1 




0. 96E25765 


C. 94693282 


-0.96776691 


1.00067712 


2 




0. 9E625659 


0. 94693466 


-0.96776960 


— 1 .OOOo77 21 


3 




0.25778800 


0. 55401616 


-1.81344349 


1 .954425 46 


4 




C. 257 7E439 


0. 55401729 


-1.61338128 


-1.954429 85 


5 




0. OCCOOOIO 


0. 30602226 


-2.406C7412 


2. £63692 97 


6 




0. CCCOOOlO 


-0. 30602324 


-2 .40740413 


-2.6637299c 


RESULTS 


CF OPTIMIZATION WITH 2 


EXTRA DATA POINTS 


PAIR 


ft 


RES. -REAL 


RES.-IMAG 


POLE-REAL FOLE-IMAG 


1 




C.97C03027 


0. 94583414 


-0.96792409 


1.00081726 


2 




C. 97CC2941 


-C. 94583770 


-0.96793108 


-1.00081754 


3 




0.25U6247 


0. 53580191 


-1.78711241 


1.95425944 


4 




0.25167543 


-0. 53579688 


-1.78713746 


-1.95923795 


5 




0. CCCCCOlO 


0. 26710278 


-2.24090416 


2.86059552 


6 




0. OCCOOOIO 


-0. 26681661 


-2 .24551510 


-2.874173 9:: 



RESULTS 


CF OFTIMIZATICN WITH 4 


EXTRA DATA PCINTS 




PAIR # 


RES. -REAL 


RES.-IMAG 


FCLE-REAL 


FOLE-IM 


AG 


1 


0.97144662 


0. 94762C10 


-0.969C4247 


1.C00805 


14 


2 


C. 97144593 


-0. 94762088 


-0.969C5146 


-1.000305 


47 


3 


0.25997577 


0. 54275341 


-1.80291116 


1.955CP2 


6c 


4 


0.25997714 


-0. 54275353 


-1.80292514 


— 1 . 955o42 


22 


5 


C. CCCOOOlO 


0. 27664033 


-2.30082557 


2 .86c760 


95 


6 


0. OCCOOOIO 


-0. 27744624 


-2.27393152 


—2 .664306 


62 
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RESULTS CF OFTiWIZATICN WITH 6 EXTR;J OAT^ PCINTS 



PAIR 


4 


RES.-REAL 


RES .-IMAG 


POLE-REAL 


1 




C.97C/4646 


0.94593786 


-0.96820496 


1 


2 




0.97C74740 


-0.94593893 


-0.96820016 


-1 


3 




0.25317420 


0.53401680 


-1.78619717 


1 


4 




0. 25317481 


-0.53401801 


-1. 78618621 


-1 


5 




o.occaoolo 


0.26390985 


-2.23469098 


2 


6 




C. OCCCOOlO 


-0.26397500 


-2.23254981 


—2 


RESULTS 


CF OPTIMIZATICN RITF 8 EXTRA DATA PCINTS 


PAIR 


4 


RES .-real 


RES.-IMAG 


POLE-REAL 


1 


1 




C. 97233004 


0. 94707234 


-0.96927693 


1 


2 




C. 97235259 


-C. 94707302 


-C.969C5353 


-1 






0.24416111 


0. 53457339 


-1.78489285 


1 


4 




C. 24415737 


-C. 53465427 


-1.79548721 


-1 


5 




O.CCCCOOlO 


0. 20134605 


- 1.70725522 


2 


6 




C. CCCuOOlO 


-0. 21428806 


-1.61287861 


-2 


RESULTS 


CF OPTi^'IZAT ICN kITF 10 EXTRA DATA PCINTS 



POLc-lM AG 



.95473002 
.954731 2i 
.86 34oo 62 
.8o 625 1 7 1 



FCLE-IHAG 

.00037271 

.00036402 

•93362845 



PAIR n RES. -REAL 

1 C. 97265329 

2 C. 97266290 

3 0.24466698 

4 C. 24466604 

5 C.OOCCOOIC 

6 C.OCCCCQIO 



RES.-IMAG 
0. 94567240 
-C. 94587709 
0. 52883260 
-0. 52886014 
C. 25592515 
-0. 15339728 



POLE-REAL 

C.97CC6813 

•0.96753425 

1.7785504? 

■1.78467362 

1.O9220053 

1.60711291 



PCLE-IMAG 
1.00100830 
1.00039816 
1 .955545 55 
1.950865 97 
2.78203752 
-3.01800673 



41 



I 




Time in nsecond 



Figure 5.1. 



Estimated E(k) for Signal 1 of 15 dB SNR. 
(32 data points from late time region) . 
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Figure 5.2. 



Estimated E(k) for Signal 1 of 15 dB SNR. 
(48 data points from late time region) . 
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Figure 5.3. Pole and Residue Plot for Signal 1 of Noise Free. 
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Figure 5.4. Pole and Residue Plot for Signal 1 of 30 dB SNR. 
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Figure 5.5. Pole and Residue Plot for Signal 
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Figure 5.6. Pole and Residue Plot for Signal 2 of Noise Free. 
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Pole and Residue Plot for Signal 2 of 30 dB SNR 
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Figure 5.8. Pole and Residue Plot for Signal 
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Figure 5.9. Pole and Residue Plot for Signal 3 of Noise Free 
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Figure 5.10. Pole and Residue Plot for Signal 3 of 30 dB SNR. 
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. 11 . Pole and Residue Plot for Signal 3 of 15 dB SNR 



C. PERFORMANCE EVALUATION 

Unlike the traditional methods, this method revealed 
itself as one that can be applicable to a general scattering 
problem, including the early time data. A potential superiori- 
ty of this method over the traditional methods is found in its 
capability of handling the early time signal that has relative- 
ly large energy content and less sensitivity to the fixed noise 
level . 

Regarding the capability of handling the early time signal, 
the performance of this method seems to be very high compared 
with the traditional methods which do not model that signal 
generically. But, two points need to be considered when we 
ensure this high degree of performance. One important point 
is that we be able to estimate properly the number of poles to 
optimally ask for. This will depend upon the SNR of the data, 
increasing generally with increasing SNR. The actual physical 
scattering data has an infinite number of poles but, because 
the energy content of natural modes tends to drop off rapidly 
beyond the bandwith of the excitation signal the higher order 
poles are ignored, blanding with the noise present. Another 
point is that the initial data window has to contain only the 
late time signal data. 

As the results of the tables in the previous pages indi- 
cate, converging characteristics of the value of parameters 
is not blocked by a high level of noise pollution. It has 
been observed that the computed value of poles are likely to 
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converge closely to the true values that are listed in Tables 
I through III either by choosing the number of poles as small 
or by increasing the size of the data window. Even though it 
is hard to prove this optimality of this method (due to the 
non-linearity of the scattering transient response signal 
model) , handling of the early time signal seems to be cast 
into a simpler problem space and that optimality may be ex- 
plained by observing the parameters at every processing step. 
Namely, the parameter estimation of this method is less 
influenced by the unreliable points than the more accurate 
ones. Figure 5.1 and 5.2 and the results in the previous 
section explain the reason why the above assertion is valid. 

We see that the estimated value of parameters through optimi- 
zation is not going to worse as we include the early time 
data points (these data points can be regarded as the unrelia- 
ble data when we extrapolate the sum of exponentials of the 
late time region into those of the early time region) into the 
data window what we see. But, as we are adding more data 
points of the late time region, improved values of parameters 
are obtained through processing, as in the results and Fig. 
5.2. Although this is true, we still see that the plot of 
E(k) in Fig. 5.2 does not provide sufficient data to confirm 
that an improvement of accuracy depends only on the number of 
late time data points, because those data points have a lower 
level of SNR locally. Now, let us get back to Tables IV, VII, 
and X. For the signal inputs whose SNRs are assumed to be 
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infinite (noise free) , we have a set of parameters whose 
values are not converging to the true values as we add an 
increased number of unknowns. These results ensure us that 
the increased computer round-off errors obtained by more 
data processing (more parameters to extract) will produce 
divergence. Even initially with low SNR, the increased 
signal energy in the additive early- time data may improve 
the results until, by adding enough extra data points and 
parameters, round-off error will eventually catch up with 
the most optimized one and beyond a specific point the 
results will begin to diverge. So that, in a general sense, 
we can say that the improvement in accuracy may be obtained 
either by the number of data points or by making use of the 
data points of high SNRs. If we consider a signal having 
high damping coefficients we can not use the signal data 
what are likely to be hidden by noise. The technique used 
here in processing the data proceeds by sliding the first 
data point along the time-axis towards the origin until the 
data window expands its range up to time-origin by adding 
the unknown parameters E (k) 's which is the same as the number 
of data points in the early time region. 

Figures 5.3 through 5.11 shows how rapidly the parameters 
we want to derive converge to the true values, within a con- 
vergence bounds. Here, all the plots were reconstructed from 
the results. 
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VI . SUMI4ARY 



A modified least-squares minimization method has been 
developed and tested against synthetically generated signal 
data using a non-linear parameter optimization algorithm. 

The effect of varying the number of data points either in 
the region of early time or in the late time under various 
typical noise environments has been studied, defining the 
criteria for using this algorithm for the extraction of poles 
and residues under relatively heavy background noise, as well 
as for the time varying residue signals modeled by Morgan 
[Ref . 4 ] . 

In an attempt to use this method for signals with high 
damping coefficients, 3 synthetically generated signals were 
used as the inputs (Appendix C) . Extraction of poles from 
the direct synthetic signal waveform was done through non- 
linear numerical evaluations, not by the Prony-based evalua- 
tions which were used in the traditional method. The 
parameter estimating algorithm for this non-linear signal 
model had worked successfully, providing users with optimized 
parameters that were converging to the true values, within 
specified convergence bounds. It was found that although we 
have the early time signal data having very high SNR, those 
data can not contribute to improving accuracy by simply 
increasing the number of extra data points and extra unknown 
parameters . 



56 



In order to get more accurate values of poles, we have 
to have an increased base of late time signal data, which 
have relatively low SNR. 

It has been shown that this method is "robust” even under 
heavy noise conditions. But three basic requirements have to 
be met for using this method as a general methodology; the 
optimal number of requested poles has to be known in advance 
to processing, reasonable initial estimation of parameters 
has to be made, and the transition time from early to late 
time signal models has to be known a priori. 
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APPENDIX A 

SYNTHETIC DATA GENERATION PROGRAM 
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APPENDIX B 

NON-LINEAR PARAMETER OPTIMIZATION PROGRAM 
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appendix c 

DATA ?;TGNAL PLQT.q 



p 




Figure C.l. 



Decomposition of Signal 1 
(Noise Free ) 
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Figure C.2. Decomposition of Signal 1 

(Noise Free) 
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Figure C.3. Decomposition of Signal 3 

(Noise Free) 
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Amplititude 




Figure C.4. Reconstruction of Signal 1 from The 

Computed Poles and Residues ( SNR=30d b ) 
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Ampiititude 




Figure C.5. Reconstruction of Signal 2 from The 

Computed Poles and Residues(SXR=30db) 
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Figure C.6. Reconstruction of Signal 3 from The 

Computed Poles and Residues ( SNR = 30d b ) 
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